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Here, we document in-vivo bite forces recorded from wild piranhas. Integrating this empirical data with 
allometry, bite simulations, and FEA, we have reconstructed the bite capabilities and potential feeding 
ecology of the extinct giant Miocene piranha, Megapiranha paranensis. An anterior bite force of 320 N from 
the black piranha, Serrasalmus rhombeus, is the strongest bite force recorded for any bony fish to date. 
Results indicate M. paranensis' bite force conservatively ranged from 1240-4749 N and reveal its novel 
dentition was capable of resisting high bite stresses and crushing vertebrate bone. Comparisons of body 
size-scaled bite forces to other apex predators reveal S. rhombeus and M. paranensis have among the most 
powerful bites estimated in carnivorous vertebrates. Our results functionally demonstrate the extraordinary 
bite of serrasalmid piranhas and provide a mechanistic rationale for their predatory dominance among past 
and present Amazonian ichthyofaunas. 

The evolution of gnathostome jaws, along with bite forces that can capture and masticate active prey is a key 
functional innovation underlying the diversification of early Devonian vertebrates 1 . As a result of their 
fundamental importance in expanding predatory niches and promoting the success of vertebrates, jaws and 
bite forces in living and extinct species have been repeatedly modeled through lever and linkage mechanics, bite 
simulations, and 3-D finite element analyses (FEA) 2 6 . However, few theoretical studies have validated their 
models with empirical measurements from living taxa 7JJ . 

In-vivo experiments in the field to elicit and record ecologically realistic biting behaviors for predatory species 
are rare, dangerous, and difficult to perform 10 . Among bony fishes, piranhas (Serrasalmidae) represent an ideal 
group of predatory vertebrates in which to investigate the evolution of extreme biting capabilities because of their 
aggressive nature, relatively small size, and accessible populations. While anecdotes of piranha-infested waters 
skeletonizing hapless victims are generally hyperbole, the effectiveness of their bite is not. Even at their small 
body sizes, diet studies indicate that piranhas will attack and bite chunks out of prey many times larger than 
themselves 1113 . 

Recent molecular evidence indicates there are three major subclades of the Serrasalmidae: 1) the carnivorous 
piranha-clade, 2) the omnivorous MyZews-clade, and 3) the herbivorous pacu-clade 14 . Within the piranha-clade, 
feeding ecology varies from the typical flesh and fin eating forms of Serrasalmus and Pygocentrus spp., to the 
highly specialized lepidophagous (i.e. scale-eating) Catoprion mento. Tooth shape and dentition patterns in 
Serrasalmidae demonstrate a strong functional relationship to diet 1113 . Carnivorous piranhas typically have jaws 
lined with a single row of 6-7 serrated, multi-cusped, triangular, blade-like teeth. Dentition and tooth morpho- 
logy have long been used to classify serrasalmid species 15,16 . In 2009, Cione et al." described a new extinct 
serrasalmid species of the Upper Miocene from the Parana geological formation in Argentina. Megapiranha 
paranensis is classified as a new genus of a giant, piranha-like species described from a single fossilized premaxilla 
jaw bone fragment that had a set of three triangular teeth set in a zig-zag pattern. The fossil teeth are also 
morphologically distinct. They have labio-lingually compressed pointed cusps with finely serrated cutting edges 
similar to a shark. In contrast, the mid tooth expands into a broad lingual shelf that is anchored to the jaw with a 
robust circular base (Fig. 1, inset). These morphological synapomorphies with extant serrasalmid species place 
M. paranensis as sister taxa to the carnivorous "piranha-clade" and a more distant intermediate relative to the 
herbivorous "pacu-clade" (SI Fig. 1). 

The dentition pattern and tooth shape of the M. paranensis fossil indicates an intermediate morphology 
capable of both slicing soft flesh and crushing hard prey. However, the diet of this giant piranha species still 
remains a mystery. The Miocene epoch is renowned for its gigantism in Neotropical aquatic flora and 
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Figure 1 | Skull anatomy of S. rhombeus (photo by SH) and fossil teeth of M. paranensis (inset, reprinted from Cione et al.). Teeth of S. rhombeus and 
M. paranensis share similarly shaped cusps that form labio-lingually compressed triangular blades. Arrow (inset) demarks fine serrations on M. paranensis 
tooth for slicing flesh. Homologous landmarks were measured as length of labial base (mm) for the 5 th premaxillary tooth (A). Allometric relationship 
describing growth of the 5 lh tooth relative to body size, measured as total length (mm) in the black piranha, S. rhombeus (B). Tooth length scales 
isometrically with increasing body size (slope = 1.02). Data were log 10 transformed and fitted to a least-squares linear regression model: Log 10 (tooth) = 
1.02Log 10 (TL) — 1.91 (r 2 = 0.96, F ( U4) = 339.8; p < 0.0001). Using this regression equation, a conservative body size of 71 cm TL was back-calculated 
for M. paranensis and was subsequently used to estimate its lower bound of bite forces (Fig. 2A). 



fauna 18,1 '. Thus, it is reasonable to assume the food resources 
available to Megapiranha would likely have required jaw forces 
and dental weaponry capable of capturing and processing very 
large prey. 

Here we report the first in-vivo bite forces recorded from wild 
specimens of the largest species of carnivorous piranha, Serrasa- 
Imus rhombeus, and describe the underlying functional morpho- 
logy of the jaws that power their bite. Using this extant species 
as an allometric surrogate, we are able to infer the feeding eco- 
logy of the giant Miocene piranha, M. paranensis, by recon- 
structing its bite forces and the mechanical capabilities of its 
unique dentition. Lastly, by removing the effects of body size, 
we are able to demonstrate that the biting abilities of these 
relatively diminutive serrasalmid fishes have biomechani- 
cally surpassed much larger iconic predators documented in the 
literature. 

Results 

In-vivo bite forces. Maximum bite force was tested in-vivo for 15 
specimens of Serrasalmus rhombeus ranging in body length from 205 



to 368 mm TL using a customized force gauge (SI Table 1). Bite force 
varied nearly five-fold from 67 N in a 0. 17 kg individual to more than 
320 N for a 1.1 kg specimen. To examine the ontogeny of bite 
performance in S. rhombeus, we plotted bite force against body size 
(Fig. 2A). In-vivo bite forces in S. rhombeus scale with significant 
positive allometry (slope = 2.30) demonstrating that jaw strength 
increases considerably faster than body length. 

Jaw functional morphology. Dissections of S. rhombeus' jaws reveal 
their powerful bite is generated by a massive adductor mandibulae 
muscle complex. It is made up of four distinct subdivisions: Al, A2 
lateral, A2 medial, A3 (Fig. 2B, C). The Al subdivision is a smaller 
fusiform muscle that originates on the ventral portion of the 
preopercle and extends dorso-rostrally to wrap around the 
coronoid process of the articular and insert onto the dorso-medial 
surface of the dentary. The two large A2 subdivisions span out and 
fill the entire suspensorium to make up more than 80% of the 
adductor mandibulae mass. Together the A2 subdivisions along 
with the medial A3 fuse into a thick rope-like tendon that inserts 
supra-distally into a deep Mecklian fossa on the medial lower jaw 
(Fig. 2C). 
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Figure 2 | Scaling of maximum in-vivo bite forces (N) plotted against body size (TL) for Serrasalmus rhombeus (A). Data were Log 10 transformed and 
fitted to least squares linear regression: Log 10 (N) = 2.30Log 10 (TL) — 3.469 (r 2 = 0.77, F (lil4 ) = 43.7; p< 0.0001). Predicted maximum bite force for a 71 
cm TL M. paranensis is also plotted (dashed arrow). Jaw functional morphology revealed a massive adductor mandibulae muscle complex (AM) that 
rotates the lower jaw postero-dorsally (white arrow) via a robust rope-like tendon (B, C). Lower jaw mechanics depict a highly modified 3 rd class lever 
where the closing mechanical advantage (Li/Lo) amplifies AM muscle force transmission by 50% to 150% from the jaw tip to the posterior teeth (C). 
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Table 1 | Comparison of in-vivo and theoretical bite forces (BF) derived from adductor mandibulae (AM) physiology and lower jaw 
morphometries for Serrasalmus rhombeus (N = 6) 



TL (mm) 


Body Mass (g) 


AM(g) 


AM CSA (cm 2 ) 


MA (jaw tip) 


In Vivo Max BF (N)* 


Anterior BF (N)* 


Mid-Jaw BF (N) 


365 


1 130 


24.7 


21.0 


0.51 


306.93 


31 1.3 


631.2 


368 


1 120 


26.5 


23.5 


0.46 


320.27 


287.6 


559.7 


295 


600 


14.5 


14.1 


0.47 


235.76 


156.9 


269.7 


322 


850 


18.9 


18.2 


0.51 


155.69 


218.4 


402.3 


210 


220 


5.7 


7.7 


0.48 


84.51 


72.3 


123.6 


240 


280 


6.3 


6.9 


0.47 


97.86 


76.2 


139.7 



'BF not significantly different: ANCOVA (t-ratio d f . 2,1 1 = 0.65, p = 0.53 1 ). 



The mechanical advantage (MA) of the lower jaw in vertebrates 
determines the proportion of adductor muscle force that is transmit- 
ted to the bite and is a function of the ratio of the in-lever length (Li = 
the distance from the jaw joint to where the adductor muscles insert 
onto the mandible) to the out-lever length (Lo = the distance from 
the jaw joint to the distal-most tooth tip). In Serrasalmus rhombeus, 
the extreme anterior insertion of the A2/A3 tendon creates a high 
mechanical advantage for even the most rostral teeth (MA ~ 0.5) 
(Fig. 2C, Table 1). For practical comparison, if the prey is bitten 
farther back along the lower jaw at the 4 lh tooth position, the shorter 
Lo increases the MA to 1.0 and approaches 1.5 at the most posterior 
teeth. 

To investigate the effects of having a massive adductor mandibulae 
complex in combination with a high MA on piranha bite force, we 
ran a 2D computer simulation that uses jaw morphometries and 
adductor muscle physiology parameters to predict bite forces in six 
individual S. rhombeus covering a range of body size. As expected 
from the jaw lever mechanics, simulation results demonstrate bite 
force nearly doubles from the anterior most teeth to the mid-jaw 
position (Table 1). For one of our largest specimens (TL = 
365 mm), predicted mid-jaw bite force (where MA = 1.0) was cal- 
culated to be as high as 631 N, more than twice the in-vivo bite force 
recorded for this individual. However, an analysis of covariance 
comparing simulation results with our empirical data revealed the 
predicted anterior bite forces derived from jaw morphometries are 
not significantly different from the maximum in-vivo bite forces 
recorded for these specimens (Table 1). 

Body size and bite force reconstruction. Estimates for Megapiranha 
paranensis were back-calculated from allometric relationships of 
tooth growth and in-vivo bite force for S. rhombeus. In contrast 
to previously published estimates (—73 kg, 128 cm TL) 17 , our 
allometric analysis of M. paranensis' fossil premaxilla suggests a 
more conservative body size: —10 kg, 71 cm TL (Fig. IB). Bite 
force predictions for M. paranensis using our calculated body size 
and the previous published estimate ranged from 1240 to 4749 N, 
respectively, at least a 4-fold increase in bite force over the largest 
S. rhombeus (Fig. 2A). 

Bite simulations and fea of fossil teeth. Bite simulations using a 
bronze-alloy metal replica of the M. paranensis jaw fossil were 
conducted to examine the potential for comminution of bony 
materials. These indentation trials tested the fossil teeth's ability to 
penetrate ecologically relevant bony prey of varying thickness. Using 
the predicted bite force range (e.g. 1245-4448 N), Megapiranha 's 
teeth penetrated the thick cortical layer (5.94 mm) of a bovine 
femur in a primarily linear fashion generating piercing inden- 
tations from 1.0-2.3 mm deep (Fig. 3). Further simulation tests on 
aquatic turtle carapace and dermal scales from armored catfishes 
repeatedly resulted in catastrophic punctures at much lower bite 
forces (range = 66.7-889.6 N; SI Movie 1). 

To better understand the functional significance of the uni- 
que shape of the M. paranensis premaxillary tooth, we used FEA 



programs to simulate the predicted loads generated during biting. 
We applied these loads to models of the M. paranensis premaxillary 
teeth, as well as to homologous teeth of two closely related serrasal- 
mid relatives that cover the spectrum of tooth shapes and feeding 
ecologies (Fig. 4). We examined both Von Mises stress, which acts as 
a predictor of the likelihood of failure, and Total Strain Energy, which 
predicts resistance to deformation and thus efficiency of the shape at 
transmitting bite forces. Using these two metrics, we examined the 
effects of different tooth shapes under the high loading regimes that 
would be experienced when biting into hard bony materials. FEA 
results indicate when force is applied at the cusp of the tooth, stresses 
are concentrated around the loading area in all three morphs. In 
addition, stresses are distributed across the body of the Piaractus 
brachypomus tooth, onto the medial face. Distribution patterns 
change when force is applied over the upper third of the tooth as if 
it were puncturing through a hard bony material. Across all tooth 
types, stresses are concentrated around the lowest point of loading. 
However, stresses are distributed through the body of the tooth in 
both the P. brachypomus and M. paranensis and are channeled to the 
base of the tooth. In contrast, stresses in the Pygocentrus nattereri 
tooth remain concentrated in a broad ring around the middle of the 
tooth. The magnitudes of Von Mises stresses and Total Strain Energy 
are consistently higher in the sharper P. nattereri and M. paranensis 
tooth morphologies compared to the blunter shape of P. brachypo- 
mus indicating a greater likelihood of structural failure and deforma- 
tion while puncturing hard materials at these forces (Fig. 4; Table 2). 

Discussion 

An in-vivo bite force of 320 N in Serrasalmus rhombeus is the stron- 
gest yet recorded for any bony or cartilaginous fish to date 7 , and is 
nearly three times greater than the bite force of an equivalent size 
American alligator 8 . Positive allometry of bite force is common 
in carnivorous and durophagous bony and cartilaginous fishes that 
excise chunks of flesh or fracture hard-shelled prey by crushing 3,20 22 . 
Our analysis of wild piranhas fits this trend, and the in-vivo data 
demonstrate that S. rhombeus can bite with a force more than 30 
times its weight, a remarkable feat yet unmatched among vertebrates. 
However, it should be noted that there was some variation in our in- 
vivo bite force data among similar sized individuals that was likely 
caused by some individuals under- performing due to stress or 
fatigue. At least one fish (TL = 310 mm) appears to be an outlier 
in the regression depressing the slope of the relationship. Removal of 
this individual increases the slope to 2.48 and improves the r-square 
of the relationship to 0.85. Thus, taking this outlier into account, the 
significant positive allometry of bite force presented here should be 
considered a conservative lower bound for maximal performance in 
S. rhombeus. 

How do black piranhas achieve these powerful bites at such rela- 
tively small sizes? The answer stems from the extraordinary size of 
the adductor mandibulae complex and the efficient transmission of 
its large contractile forces through a highly modified jaw closing 
lever. Indeed, the mass of the adductor mandibulae complex makes 
up over 2% of S. rhombeus' total body mass. The highest percentage 



SCIENTIFIC REPORTS | 2 : 1009 | DOI: 1 0. 1 038/srep01 009 



4 



5000 



4500 



4000 



3500 



3000 



g 2500 



2000 



1500 





1 cm 



1000 



500 



0.00 0.50 1.00 1.50 2.00 

Indentation (mm) 




2.50 



3.00 



Figure 3 | Bite simulation curves using a metal-alloy replica of M. paranensis fossil upper jaw (left inset) indenting a bovine femur (A), aquatic turtle 
carapace (O, Phrynops sp.), and dermal scales from an Amazonian armored catfish (+, Pterygoplichthys multiradiatus). Black arrows indicate the 
predicted anterior bite force for small and large body size estimates of M. paranensis. The 3 rd and 5 th teeth punctured the cortical bone layer of the bovine 
femur from 1.0-2.3 mm (inset upper right corner). Results from the turtle carapace and catfish dermal scales reveal Megapiranha' s dentition pierces 
through these bony prey items at relatively low forces. 



yet recorded for bony fishes 23 . This large muscle mass combined with 
a short overall muscle length results in a huge physiological cross 
sectional area that is directly proportional to the amount of force it 
can produce (Fig. 2B, Table 1). Furthermore, the extreme anterior 
insertion of the A2/A3 tendon gives S. rhombeus an anterior bite 
force that is powered by one of the highest jaw closing mechanical 
advantages ever identified in fishes (Fig. 2C, Table l) 24 . Thus, unlike 
many actinopterygiian fishes where the closing MA of the lower jaw 
averages 0.27 (i.e. less than 30% of the adductor's muscle force is 
transmitted to the bite), S. rhombeus lower jaw lever can harness 50% 
of the adductor complex's force generating capacity at the tip of its 
jaws. Even more impressive, if the prey is bitten on the lateral side of 
the jaws near the posterior teeth, this musculoskeletal arrangement 
further increases the jaw's closing MA to amplify the adductor com- 
plex's torque during biting by as much as 150% (Fig 2C). Our 
simulations results demonstrate how this underlying functional 
morphology amplifies the adductor muscles' force so that even the 
smallest specimen (0.22 kg) can potentially generate a range of trau- 
matic bite forces across the jaw (e.g. 72-123 N; Table 1). The fact that 
our theoretical estimates for anterior bite force are statistically indis- 
tinguishable from the bite performances that were recorded near the 
tip of the jaws indicates the 2D model's simulation of lever mechanics 
combined with adductor muscle physiology accurately explains bite 
force capability in S. rhombeus. So, when one considers their unique 



jaw functional morphology combined with their aggressive biting 
behavior, it should come as no surprise that black piranha whether 
large or small can rapidly and efficiently excise large chunks out of 
their prey. 

Although our body size estimate of M. paranensis is considerably 
smaller than previous published accounts 17 , it is still more than triple 
the maximum body size of extant S. rhombeus 15 . It is also important 
to note that both of these estimates fall within the range of sizes for 
the largest living serrasalmid, the frugivorous Tambaqui, Colossoma 
macropomum 15 . Paleontological reconstructions using fossil evid- 
ence and inferences drawn from related extant taxa are commonly 
used to generate hypotheses about ecology, behavior, and biomecha- 
nics of extinct species 2,4,5,25,26 . Because of their close phylogenetic 
relationship, we assumed homologous jaw biomechanics, and used 
the ontogeny of S. rhombeus as a proxy to predict Megapiranha s bite 
forces. For comparison, our analysis predicts Megapiranha s bite was 
equivalent to the anterior bite force of a great white shark weighing 
over 400 kg (SI Table 3). Again, these bite force predictions likely 
underestimate Megapiranha' s maximum bite force since they are 
based on the allometry of the variable in-vivo data. Furthermore, 
the in-vivo allometric relationship represents anterior bite forces 
and does not take into account the potential for doubling bite forces 
along the lower jaw. Thus, if M. paranensis did have similar muscu- 
loskeletal architecture in the lower jaws as in S. rhombeus (Fig. 2C) 
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Figure 4 | Finite Element Analysis of serrasalmid dentition. Maximum and minimum Von Mises stresses reported for maximum predicted bite force 
(scaled relative to predicted bite forces in Megapiranha); distribution of stresses remains the same regardless of magnitude of applied force. Warmer colors 
indicate higher concentrations of Von Mises stress and cooler colors reflect lower stress concentrations. Higher concentrations of Von Mises stress 
indicates a greater likelihood of failure in that region of the tooth. 



but scaled up to its larger body sizes, then the mechanical advantage 
at its mid-jaw (MA = 1.0) could have produced bite forces as high 
as 2480 N to 9498 N. With these amplified estimates, a 73 kg 
Megapiranha s biting attack would have had the same ferocity of a 
3000 kg great white shark (SI Table 3f. 

What Megapiranha fed on with these powerful bites is uncertain. 
Tooth morphology commonly reflects feeding ability in fishes, with a 
strong functional relationship to diet in the Serrasalmidae 11,16,23 . 
However, the unique dentition of M. paranensis presents a paradox, 
exhibiting traits from relatives of both the herbivorous pacu-clade 
and the carnivorous piranha-clade 14,17 . Megapiranha teeth have 
robust circular bases with a broad lingual shelf that tapers into short 
serrated triangular blades at the crown, making them potentially 
suited for processing both soft and hard prey materials (Fig 2, inset). 
Did Megapiranha use their strong jaws and novel dentition to not 



only slice into their prey's compliant flesh, but also crush through 
their stiff bones? Our bite simulations indicate Megapiranha' s hybrid 
teeth were indeed capable of transmitting sufficient bite pressure to 
generate micro-fractures in robust cortical mammalian bone as well 
as complete mechanical failure of thinner vertebrate bony tissues 
(Fig 3). For comparison, even modern carnivorous piranhas species 
with less robust tooth shapes are capable of some bone comminution 
as they often feed by taking bites out of the bony fins of fishes 1 '~ 13 . In 
fact, there are also documented cases of S. rhombeus biting off and 
consuming human phalanges 27 . Thus, the combination of applying 
the predicted bite forces from the allometry of S. rhombeus with the 
fossil replica in our penetration tests demonstrates the tooth shape of 
M. paranensis would have facilitated an osteophagous diet of turtles 
and armored catfish and to a lesser extent limb structures of larger 
terrestrial mammals during the Miocene Epoch. 



Table 2 Finite Element Analysis of applied forces (N) 


and output for serrasalmid teeth 












Maximum Von Mises 


Minimum Von Mises 




Relative Total Strain 


Species 


Type of Load 


Stress- Scaled (N) 


Stress (Pa) 


Stress (Pa) 


Strain-Scaled (N) 


Energy (J) 




Cusp tip 


677.52 


1.18E+1 1 


4.79E-04 


183.333 


2.64E-03 


Pygocentrus nattereri 




190.937 


3.32E+10 


1.35E-04 


51.667 


2.1E-04 


Distal third 


677.52 


6.21E+09 


1.82E-04 


183.333 


5.94E-04 






190.937 


1 .75E+09 


5.12E-05 


51.667 


1.7E-03 




Cusp tip 


4400* 


2.31E+1 1 


4.06E-03 


4400* 


2.85E-01 


Megapiranha paranensis 




1240* 


6.52E+10 


1.14E-03 


1240* 


2.26E-02 




Distal Third 


4400* 


6.76E+09 


5.75E-04 


4400* 


5.91E-02 






1240* 


1.91E+09 


1.62E-04 


1240* 


4.7E-03 




Cusp tip 


476.59 


3.79E+10 


5.27E-03 


305.556 


1.86E-03 


Piaractus brachypomus 




134.312 


1.07E+10 


1.48E-03 


86.1 1 1 


1.47E-04 


Distal Third 


476.59 


2.14E+09 


5.59E-03 


305.556 


3.2E-05 






134.312 


6.03E+08 


1.57E-03 


86.1 11 


4.03E-04 


*No scaling. 
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From the FEA models, the distribution of Von Mises stresses in the 
M. paranensis tooth revealed hybrid patterns that reflect those 
observed in both the sharp triangular blades of the carnivorous red 
belly piranha, Pygocentrus nattereri, and the blunt cuspids of the 
durophagous pacu, Piar actus brachypomus (Fig. 4). When loaded 
at the tip, both M. paranensis and P. nattereri exhibit Von Mises 
stress distributions similar to FEA models of shark teeth puncturing 
soft prey 28 . In contrast, once the tip has presumably pierced into the 
bony material and compressive loading regimes are focused on the 
upper third of the tooth, the shape of M. paranensis distributes Von 
Mises stresses like the half-domed mollariform shape of P. brachy- 
pomus. The pattern reveals a reduction of Von Mises stresses in the 
broader sloped middle portion of the cusp where increased enameled 
surface area would come in contact with the bony prey material 
material after initial penetration (Fig. 4). The distribution patterns 
of Total Strain Energy (not shown) also closely mirrored the Von 
Mises stresses, indicating higher efficiency and lower deformation in 
the same region. Contrary to our initial hypothesis, the hybrid mor- 
phology of the M. paranensis tooth revealed the highest Von Mises 
stresses and Total Strain Energies of the three species, though the 
differences are slight (Table 2). This suggests that there are perform- 
ance tradeoffs when combining shapes for piercing and crushing in 
the same tooth. However, while relative stresses and strain energies 
are high in M. paranensis, their distribution patterns reflect tensile 
strains, not only similar to those we found in P. brachypomus, but 
also seen in FEA modeling of mammalian cingulum 29 . A partial 
cingulum in mammals creates a morphologically similar shelf that 
reduces tensile strains patterns in the tooth by as much as 21% after 
penetrating the prey. This enhanced strain resistance is apparently 
achieved by increasing the total surface area of the shelf that comes in 
contact with load. The M. paranensis tooth is apparently able to 



overcome these scaled shape limitations in a similar way by being 
absolutely larger than the teeth of either P. nattereri and P. brachy- 
pomus and expanding the surface area of the lingual slope where 
crushing forces would be most efficiently transmitted. Thus, the 
novel shape of Megapiranhas dentition appears functionally adapted 
to not only focus stress at the tip while puncturing flesh like the teeth 
of a red belly piranha, but also effectively spread impact stresses 
through the thicker broader base for cracking through bony materi- 
als similar to a nut- crushing pacu. 

From these integrated simulations and FEA models, we surmise 
feeding on significantly hard prey items was at least biomechanically 
possible for M. paranensis. Confirmation of true osteophagous beha- 
vior on large vertebrate prey, however, still requires the discovery of 
Miocene fossils with bite marks that can be attributed to its dentition. 
Even so, our results indicate that the predicted bite pressures and 
hybrid tooth morphology of Megapiranha would have granted access 
to a vast menu of large prey resources present during the Miocene 18 . 

The reconstructed bite forces of M. paranensis are impressive, but 
how do they compare to other vertebrate mega-predators? Using 
fossil evidence and indentation simulations comparable to those 
presented here, the bite force of Tyrannosaurus rex was estimated 
to be over 13,400 N 26 , almost three times the largest M. paranensis 
estimate. However, much of this difference is due to scaling effects 
attributable to T. rex's over 100-fold greater body mass. To control 
for the effects of body size, we used published theoretical estimates 
for anterior bite force, and calculated normalized bite force quotients 
(BFQ's) 25 for nine living and extinct bony and cartilaginous fishes 
(Fig. 5; SI Table 3). These species vary widely in their phylogenetic 
and geological origins, but all are apex predators that feed on large 
prey with similar feeding strategies and powerful jaw mechanics. 
Correcting for body size differences reveals that both the living 
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Figure 5 | Bite Force Quotients of living and extinct apex predatory fishes. Comparison of anterior bite forces among apex fish predators using 
calculated Bite Force Quotients (BFQ). BFQ's correct for absolute differences in bite forces attributable to drastic differences in body size by comparing 
mass specific residuals. Species with generally average bite strengths for their body size have a BFQ of 100. Black and white bars represent the lowest and 
highest quotients from small and large body size estimates, respectively. Barracuda and shark images © Gray Taxidermy. Dunkleosteus terrelii image © 
Karen Carr. 
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Serrasalmus rhombeus and the extinct Megapiranha paranensis have 
among the most powerful bites in carnivorous fishes, living or 
extinct. Indeed, for its relatively diminutive size, Megapiranha para- 
nensis' bite dwarfs other extinct mega-predators, including the 
enormous whale eating Carcharodon megalodon and the monstrous 
Devonian placoderm, Dunkleosteus terrelli. 

The extreme biting abilities and jaw functional morphology of 
wild black piranhas presented here provide the first quantitative data 
that illustrate the effectiveness of their feeding mechanism. In addi- 
tion, if our fossil reconstructions and simulations are correct, then 
Megapiranha paranensis was indeed a ferocious bone-crushing 
mega-predator of the Miocene epoch. Taken together with our body 
size corrected bite force estimates, our results for living and extinct 
species validate the fearsome predatory reputations of piranhas and 
establish this group of fishes as a pinnacle of performance in gnathos- 
tome jaw evolution. 

Methods 

Specimen collection. In August 2010, we collected live black piranha, Serrasalmus 
rhombeus (Linnaeus, 1766), while fishing along the Xingu and Iriri tributaries of the 
Amazon River basin near Altamira, Brazil. A total of fifteen fish ranging in body size 
from 205-368 mm TL were collected by hook-and-line and gill nets. To ensure 
minimal damage to the jaws and record maximal performance, we used barbless 
hooks and carried out bite force experiments as soon as possible after capture to 
reduce stress and fatigue. After the expedition, collected specimens were donated to 
the INPA Ichthyological Collection at Manaus, Brazil. 

Bite performance. In-vivo bite forces were recorded using a customized portable and 
waterproof force gauge constructed of a LCGD -250 low profile miniature load cell 
(range — 0-250 lbs) and a DP7600 digital high-speed load/strain meter (500 Hz 
sampling rate) (www.omega.com). This species readily performed multiple defensive 
bites with the transducer placed between the tips of the jaws. After each bite, peak 
force was recorded and the force gauge was removed from the jaws and tared before 
the next attempt. After bite force trials, fish weight (g) and total length (mm) were 
measured, and the fish was either euthanized in an ice bath for later anatomical 
dissections or released alive. All bite performance experiments with live wild caught 
piranhas comply with the guidelines of the Western Kentucky University 
Institutional Animal Care and Use Committee. 

Jaw functional morphology. In the field, individual 5. rhombeus were sacrificed and 
dissected to investigate the morphology of the lower jaw mechanism. For each fish, 
the main subdivisions of the adductor mandibulae (Al, A2, and A3), the major jaw 
closing muscles, were excised and weighed with an Ohaus HH-320 digital scale to the 
nearest 0.1 g. Adductor mandibulae subdivisions were carefully investigated for their 
origins, insertions, and lines of action. Musculoskeletal anatomy follows the 
nomenclature conventions of Machado-Allison 16 . For six individual fish covering a 
range of body sizes (210-368 mm TL), digital photographs of the dissected jaw 
anatomy were taken from the fish's right side with a scale bar in frame. Photographs 
were analyzed using ImageJ64 (http://imagej.nih.gov/ij) to determine landmark 
coordinates of adductor muscle insertions and origins and lower jaw dimensions. 
Morphometric data were then analyzed with MandibLever v3.5 30 . MandibLever 
calculates static and dynamic bite forces from jaw lever mechanics, non-linear length- 
tension muscle contractile properties using the Hill equation, and estimates of 
adductor muscle physiological cross-sectional area. ANCOVA was then used to 
compare the allometries of simulation results and in-vivo bite performance with body 
size as the covariate. A non- significant interaction term indicated parallel slopes 
between these two relationships. Subsequently, the ANCOVA model was rerun 
without the interaction to examine if a significant difference existed between the 
predicted anterior bite forces generated by MandibLever and the empirically 
measured in-vivo bite forces recorded in the wild. 

Scaling analyses. We summarized the mean and standard deviations of bite 
performance for each individual (SI Table 1). The peak force (N) recorded for each 
individual was used as a conservative estimate of the maximum bite performance for 
that individual. To scale bite forces in S. rhombeus, maximum bite force was regressed 
against total length (TL). The null hypothesis was isometry (i.e. slope = 2.0) in 
maximum bite force relative to body size. A second scaling relationship was also 
developed to examine the growth rate of the 5 th premaxillary tooth in S. rhombeus. 
The null hypothesis for isometric growth in tooth length was an expected slope of 1.0. 
The regression equation for this allometric relationship was then used as a proxy to 
back- calculate an estimated body size for M. paranensis using the homologous length 
for its 5 th premaxillary fossil tooth. Allometric relationships were developed using 
standard least square regression models. All data were Log 10 transformed and 
analyzed using JMP v5 by SAS Institute (2002) statistical software program. 

Reconstructing the bite of megapiranha paranensis. We have generated 
hypothetical estimates of the body size and maximum bite forces of M. paranensis by 
extrapolating predictions from allometric relationships of its close living relative, 



Serrasalmus rhombeus. S. rhombeus is the largest carnivorous species of piranha in the 
Amazon River basin growing to over 3 kg. From the morphological phylogeny, the 
monotypic M. paranensis is nested between Catoprion and the higher monophyletic 
clade of Pygopristis, Serrasalmus, and Pygocentrus 17 (SI Fig. 1). Thus, we chose S. 
rhombeus as a suitable surrogate piranha morphotype with which to reconstruct M. 
paranensis bite forces based on their grossly similar carnivorous tooth morphologies 
and close phylogenetic relationship. To reconstruct hypothetical bite forces in M. 
paranensis, we used two estimates of body size: the smaller 710 mm TL, based on the 
scaling analysis of tooth growth rates in 5. rhombeus in this study (Fig. 2), and the 
larger 1280 mm TL estimate 17 . From these two body sizes, we calculated lower and 
upper bounds for maximum bite forces using the regression equation generated from 
the scaling relationship of 5. rhombeus' in-vivo bite forces (Fig. 1A). 

Bite simulations. To investigate the effects of tooth shape on fracturing vertebrate 
bone, bite simulations were conducted using a metal replica of the M. paranensis fossil 
teeth and premaxilla. The full-scale dental-grade model was forged from bronze-alloy 
by Master Craftsman Studios. Metal alloy replicas of fossil teeth provide sufficient 
yield strength and hardness to perform as an enameled tooth analog for testing bite 
indentations 4,26 . Simulation trials penetrated ecologically relevant bony prey materials 
of varying thickness (e.g. cortical layer of freshly thawed bovine femur, turtle 
carapace, and dermal fish scales). Indentation trials examined normal loading of the 
Megapiranha replica puncturing into the material. The bony materials were secured 
in a United Testing Systems SSTM-10 KN universal loading frame and indented with 
the replica at a rate of 1 mms -1 with a displacement sampling rate of 20 Hz. Each trial 
tested the range of predicted bite forces from this study (e.g. 1240-4749 N) (Fig 3). 

FEA of serrasalmid tooth shape. For M. paranensis fossil tooth comparisons to 
extant serrasalmid species, homologous premaxillary teeth were isolated from 
preserved specimens from the Field Museum of Natural History Fish Collection for 
Pygocentrus nattereri (FMNH 108184) and Piaractus brachypomus (FMNH 110198). 
Three dimensional tooth models were generated for each species using a Next Engine 
360" Digital Laser Scanner. Individual teeth were then isolated from laser scans using 
MeshLab (vl.3.1 Visual Computing Lab - ISTI - CNR), which was also used to fill 
holes and simplify the meshes. Tooth meshes were then uploaded into Amira (v5.2.2 
Visage Imaging), where they were converted from shells to solids. Solid models were 
imported into Marc Mentat (2010.1.0 MSC Software Corporation) for analysis, and 
elements were set to be solid tetras. There has been little work done on non- 
mammalian tooth materials properties, so we used an average of the Young's Moduli 
found for shark species (7.05 e 08 N/m 2 ) 31 and a Poisson's ratio of 0.3, a standard for 
mineralized materials. Nodes at the base of each tooth model, where the tooth would 
articulate with the underlying bone, were held motionless while the rest of the model 
was allowed to move and rotate without constraint. Each model was subjected to two 
different types of loading at both the maximum and minimum predicted bite force for 
Megapiranha paranensis. To compare between teeth of different sizes, we scaled 
applied forces such that the ratio offeree to approximate surface area would be similar 
between treatments 32 . For the first treatment, to mimic the point at which the tooth 
first interacts with a hard surface, we loaded a small area at the very tip of the tooth 
cusp. In the second treatment, to simulate the forces exerted when the tooth has been 
driven into the prey, we loaded each tooth to a depth approximately equal to one third 
of the tooth height. This depth was based on penetration distance taken during the 
indentation trials. We measured the Von Mises stresses distributed in each model, as a 
proxy for structural failure; higher Von Mises values means a greater likelihood of 
structural failure. Stress distribution patterns in each tooth were also compared across 
the different treatments. To measure Total Strain Energy, tooth models were 
imported into COMSOL Multiphysics (version 4.3.0.151) and subjected to the same 
load treatments. The force used for each tooth was scaled so that the ratio of applied 
force to approximate volume was equivalent for comparison between models 32 . 

Mega-predator comparisons. Bite Force Quotients (BFQ) were calculated following 
the method of Wroe et al. 25 . Briefly, maximum anterior bite forces (N) and body sizes 
(g) for each of seven additional predatory species were gleaned from the literature (SI 
Table 2). Where possible, small and large bite force estimates for each species were 
included along with the low and high range of bite force measurements for S. 
rhombeus and estimates for M. paranensis from this study for a total of sixteen data 
points. Data were then log 10 transformed and all species bite forces were regressed 
against body size to generate a Mega-predator linear regression equation (Log 10 
(Predicted BF) = 0.59 Log 10 (g) + 0.16; r 2 - 0.84, P < 0.0001). From this Mega- 
predator linear regression 'predicted' bite forces for each species body size were 
calculated, and the size corrected BFQ for each species was generated with the 
following equation: (Anterior Bite Force/Predicted Bite Force) * 100. 
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